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The critical behavior of the model with iV-vector complex order parameter and 
three quartic coupling constants which describes phase transitions in unconventional 
superconductors, helical magnets, stacked triangular antiferromagnets, superfluid 
helium-3 and zero-temperature transitions in fully frustrated J osephson-j unction 
arrays is studied within the field-theoretical renormalization group (RG) approach 
in three dimensions. To obtain qualitatively and quantitatively correct results per- 
turbative expansions for /3-functions and critical exponents are calculated up to 
three-loop order and resummed by means of the generalized Pade-Borel procedure. 
Fixed point coordinates, critical exponents values, RG flows etc. are found for phys- 
ically interesting cases N = 2 and N = 3. Critical (marginal) values of iV at which 
the topology of the flow diagram changes are determined as well. It is argued, on 
the base of several independent criterions, that the accuracy of the numerical results 
obtained is about 0.01, an order of magnitude better than that given by resummed 
two-loop RG expansions. 

In most cases the systems mentioned are shown to undergo fluctuation-driven 
first-order phase transitions. Continuous transitions are allowed in hexagonal d— 
wave superconductors, in planar helical magnets (into sinusoidal linearly-polarized 
phase), and in triangular antiferromagnets (into simple unfrustrated ordered states) 



1 



with critical exponents 7 = 1.336, v = 0.677, a = -0.030, = 0.347, r\ = 0.026 
which are hardly believed to be experimentally distinguishable from those of the 
3D XY model. The chiral fixed point of RG equations is found to really exist and 
possess some domain of attraction provided N > 4. So, magnets with Heisenberg 
(N = 3) and XY-like (N = 2) spins would not demonstrate chiral critical behavior 
with unusual values of critical exponents; they can approach the chiral state only via 
first-order phase transitions. 
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I. INTRODUCTION 



The renormalization group (RG) approach in three dimensions proved to be very efficient 
when used to study the critical behavior of simple 0(n)-symmetric models. Calculations 
of higher-order RG expansions for field-theoretical /^-functions and critical exponents com- 
bined with proper resummation of the series obtained resulted in the estimates of critical 



This approach enabled one to give quantitatively correct description of critical behavior of 
more complex systems possessing two quartic coupling constants in their Landau-Wilson 
Hamiltonians. We mean the impure Ising model |3|-[5|], the cubic model @||, and the mn- 
model |7| . Moreover, the method turned out to be powerful enough even in two dimensions 
as was shown by comparison of the approximate results obtained on the base of four-loop RG 
expansions resummed by means of Pade-Borel-like technique with their exact counterparts 
known for (exactly solvable) 2D Ising and impure Ising models Jl],@J|. 

On the other hand, there are numerous models with more than two quartic coupling 
constants which describe phase transitions in a variety of systems. Such models, however, 
being extensively studied in the frame of RG approach were actually treated only within the 
lowest-order (one- and two-loop) approximations which are known to lead to rather crude 
quantitative and, sometimes, to contradictory qualitative results. 

The aim of this paper is to study the static critical behavior of the three-dimensional 
model with three quartic coupling constants on the base of three-loop RG series resummed 
in the way which provides proper physical predictions and accurate numerical estimates. 
As far as we know, this is the first attempt to get reliable, numerically correct results for a 
complicated 3D field-theoretical model from higher-order RG expansions || . The Landau- 
Wilson Hamiltonian of the model is as follows: 



where ip a is a complex vector order parameter field, a, (5 = 1, 2, . . . , N, a bare mass squared 
rrig being proportional to the deviation from the mean-field transition point (line). 

This model describes critical phenomena in a plenty of substances. Their list includes 



exponent values which nowadays are referred to as most accurate (canonical) numbers |T],H . 
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tetragonal, hexagonal and cubic superconductors with d— or p-wave pairing || as well as 
superconductors with two — s and d — order parameters [ 10| , fully frustrated Josephson- 



junction arrays (FFJJA) at zero temperature ||11|| , stacked triangular antiferromagnets (STA) 
|Ti2| , [TB[ | , helical magnets (HM) and magnets with sinusoidal spin structures JT^|-p^[|, the A 
phase of superfluid 3 He P!7|JT^1 . The model Eq. is related also to the critical thermo- 
dynamics of type-II superconductors with short coherence length near the upper critical 
magnetic field [[II] . 



All the systems mentioned were extensively studied during the last decade and rich 
theoretical information about their critical behavior has been obtained both analytically 
and numerically. Unfortunately, the major part of these data turns out to be contradictory 
or inconclusive. To illustrate this point we overview, in brief, what was predicted for FFJJA, 
STA and HM by different people within different approaches. 

A superconductor-insulator transition in FFJJA at zero temperature produced by com- 
petition of Josephson and charging effects in the presence of quantum fluctuations is de- 
scribed by three-dimensional model Eq. (|l.l| ) with N = 2 and Vq = or wq = 0. Starting 
from 4 — e dimensions, such a transition was shown to be, within the lowest order in e, 



discontinuous [20], while (2 + e)-expansion did not enable one to resolve whether it should 
be first-order or continuous [^UJ. On the other hand, this transition was referred to as 
second-order one on the base of analysis valid to the leading order in 1/N . 

For N = 2 and N = 3, the Hamiltonian ( |1 . 1|) governs the critical behavior of STA such as 
VC1 2 , VBr 2 , CsMnBr 3 , CsVCl 3 , and of HM (Ho, Dy, Tb, (3-Mn0 2 , MnAu 2 ). In the case of 
Heisenberg spins, RG calculations in 4 — e dimensions and 1 /TV-expansion result in a first- 



order phase transition |1J-|16|,^TJ , although e-expansion predictions were believed also as 



favoring continuous one |0|. Monte-Carlo simulations in 3D seem to provide an evidence 



of continuous phase transition [12] ) |22| . RG analysis of corresponding (2 + e)-dimensional 
model proposes that systems mentioned should undergo, in three dimensions, a first-order 
transition or a second-order one with either 0(4) (not 0(6)) critical or tricritical mean-field 
exponents p3f . 

Obviously, the situation needs to be cleared up. Since the problem does not allow exact 
solution, in order to obtain reliable theoretical information one has to employ approximate 
methods with controlled or, at least, known level of accuracy. Regular RG perturbation 



theory in 3D subject to the application of Pade-Borel-like resummation technique will be 
shown to play a role of such a method. 

The paper is organized as follows. In Sec. the Hamiltonians describing the systems 
mentioned above are considered and related to the Hamiltonian ( II . 1| ) . In Sec. |TJ the renor- 
malization scheme is formulated and three-loop RG expansions for /3-functions and critical 
exponents are presented. Various resummation techniques based on the Borel transforma- 
tion and applicable to divergent (asymptotic) power series of several independent variables 
are considered and criterions for the choise of the best one are established. The specific 



symmetries of the model ( |1 . 1| ) with N = 2 are discussed in detail. They relate coordinates 
of different fixed points of RG equations to each other being a sensitive indicator of the 
quality of the approximation employed. All the numerical results obtained are presented in 
Sec. [TV} coordinates of the fixed points, critical exponent values, critical (marginal) order 
parameter dimensionalities N c at which the topology of flow diagrams changes, etc. RG flows 
are also shown in the planes where stable, within these planes, fixed points exist. In Sec. [V] 
the results obtained are applied to superconducting, superfluid and magnetically-order sys- 
tems and certain theoretical predictions are made. Particular attention is paid to what is 
known as chiral critical behavior and its relevance to real HM and STA with Heisenberg or 
Ay-like spins. Sec. |VT] containes a summary of the results obtained. Details of Pade-Borel 
resummation procedure are described in Appendix. 



II. RELEVANT SUBSTANCES AND STRUCTURES 

In this section we discuss physical systems undergoing phase transitions which are de- 
scribed by the Hamiltonian ( |1 . 1| ) . 



A. Unconventional superconductors 

These materials should be mentioned first since Eq. ( |1 . 1|) is actually an obvious general- 
ization of appropriate Ginzburg-Landau form (see, e.g. Ref. ||) with ip a being a supercon- 
ducting order parameter. For N = 2 the Hamiltonian under consideration governs a static 
critical behavior of tetragonal and hexagonal (v = 0) superconductors with <i-wave pairing, 
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while the case N = 3 corresponds to cubic p-wave materials. Heavy-fermion compounds 
such as UPt 3 , UBei3, and others are thought to belong to this class of superconductors 



24,25]. Phase transitions in thorium-doped UBei3 are well described by the phenomeno- 



logical model with two coexisting, s- and d-wave, order parameters [TT1|,|2T)|] which, in some 



limit, is reduced to that given by Eq. (LI). Moreover, since there are numerous experi- 
mental facts [0)0 anc ^ theoretical predictions [^,0] favoring non-trivial pairing modes in 
high-T c superconductors, the Hamiltonian ( |1 . 1| ) may be also relevant to the critical behavior 
of these new materials ||31|| . 

It is worthy to note that the width of critical region is large enough in high-T c super- 
conductors (see Refs. P2|JI3|| for an overview and numerical estimates) and superconducting 
fluctuations proved to be clearly seen in their thermodynamics near T c |34]-|37|. Wide fluc- 
tuation regions are also expected to exist in heavy-fermion compounds P8| , P§] . That is why 
the critical behavior of the model (|1 . 1| ) is extensively studied within the context of supercon- 
ductivity |4"0|-j4"2"|| . On the other hand, the Hamiltonian ( |1 . 1| ) has only one, isotropic gradient 
invariant, i.e. it ignores a crystallographic anisotropy of the order parameter correlation 
function which may play an appreciable role in the critical region. So, the applicability of 



the Eq. (1J.) to real unconventional superconductors is somewhat limited. The influence 
of anisotropic gradient terms on thermodynamics of these materials in the region of weak, 
Gaussian fluctuation was studied in Ref. p3| . 



B. Fully frustrated Josephson junction arrays 

Main features of JJA behavior are known to be described by the following Hamiltonian 

H = ~Y e(Jt) -EjJ^cosie.-ej-Aj) , (2.1) 

where 0j is a phase of superconducting order parameter in z-th island, 

2 j 

A tj = ^ J Adl , (2.2) 

i 

A being a vector potential of external magnetic field, and $o is a quantum of flux. Here E c 
plays a role of charging energy which is responsible for the Coulomb blockade and quantum 
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dynamics while the Josephson coupling Ej favors establishing of the global phase coher- 
ence and overall superconductivity in the system. At zero temperature superconductor-to- 
insulator transition occurs when the ratio E c /Ej exceeds a critical value. Since quantum 
fluctuations are essential in the case considered, the effective dimensionality of the system 



turns out to be equal to three: D = 2 + 1 (see, e.g. Ref. |l ![). 

If the external magnetic field B is uniform JJA is regulary frustrated with the frustration 
parameter / = (Bao) / $o> a o being an area of a plaquette. We shall consider JJA with square 
and triangular lattices in a magnetic field corresponding to / = | which are usually referred 
to as fully frustrated ones. To study their critical behavior a proper Hubbard-Stratonovich 
transformation [0,^] may be applied to the model ( |2.1| ) resulting in the Landau- Wilson 



Hamiltonian with quartic terms which are, in notations of Ref. [ II], as follows: 



m^i| 2 + ^2| 2 ) 2 -^i|^i| 2 |^2| 2 + ^Re(^2) 2 , (2.3) 

where ipi and ip 2 are complex scalar fields. In the case of square lattice u > 0, V\ — v 2 > 0, 
while for the triangular FFJJA u > 0, V\ < and t> 2 = 0. It is easy to see that Eq. ([2.3|) is 



actually identical to the quartic part of the Hamiltonian (171) for iV = 2 provided coupling 
constants are related to those standing in Eq. (|1.1| ) by 

u = u + v + w , vx = 2(v + w ) , v 2 = 2w . (2.4) 

Domains v = 0, w > and v < 0, w = correspond to the square and triangular FFJJA 
respectively. The Hamiltonian ( |1 . 1|) governs also the critical behavior of triangular JJA with 
f — l since it is known to belong to the same universality class as FFJJA with square lattice 



C. Stacked triangular antiferromagnets 

Triangular antiferromagnets which we shall deal with possess lattices consisting of trian- 
gular antiferromagnetic layers stacked in register along the orthogonal axis. In the ground 
state the spin arrangement may be thought as formed by three ferromagnetic sublattices 



with 120° angles between neighboring, within the layer, spins (see Refs. |T3|J46[1 for detail). 
The microscopic Hamiltonian modelling STA reads 



H=-JJ2 - J' E S i S i > J < • ( 2 - 5 ) 

<ij> <ij>' 

The first sum represents antiferromagnetic interactions within triangular layers which give 
rise to frustration. The second one describes interlayer coupling, the sign of J' being 
unimportant since there is no frustration along the orthogonal direction. The Hubbard- 
Stratonovich transformation followed by the expanding around the instability points and 
other standart procedures leads to effective Hamiltonian containing 

Mfc (a 2 + b 2 ) + i; fe [(ab) 2 -a 2 b 2 ] (2.6) 



as an interaction term, a, b being real n-component vector fields [13] . If then one put 



Uk = u + w , v k = Aw Q (2.7) 



Eq. (|2.6|) will immediately turn into the quartic part of the Hamiltonian (|1 . 1| ) with t> = 
and <f a = a a + ib a . The frustration may be shown to be relevant only for wq > 0; the 
opposite case, Wq < 0, corresponds to simple ferromagnetic or antiferromagnetic ordering 

m. 



D. Helical magnets 

In these magnets spins aling ferromagnetically in a plane and form spirals along the 
orthogonal axis. Such an ordering may be described by the microscopic Hamiltonian (|2.5|) 



provided the first and second sums are defined in a new manner [f[3 |. Namely, let the first 
sum represents nearest-neighbor ferromagnetic interactions, J > 0, while the second term in 
( p.5| ) describes antiferromagnetic, J' < 0, next-nearest-neighbor interactions acting along 
only one crystallographic axis. Then for ratios | J'\j J exceeding a critical value spins will be 
helically arranged along this axis. All the machinery mentioned above gives in this case just 
the same Landau-Wilson Hamiltonian as for STA. The helical ordering, however, is realized 
only if Vk > (w > 0) fL3| . For v k < & sinusoidal (linearly-polarized) spin density wave 



should occur 16 . 
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E. Superfluid helium —3 



In liquid 3 He fermionic excitations are known to form, below T c , Cooper-like pairs with 
S = L = 1. Since magnetic dipole-dipole interaction couples orbital and spin angular 
momenta to each other the superfluid order parameter possesses a symmetry 0(3) x U(l). 
This is precisely the symmetry underlying the Hamiltonian ( |1 . 1|) with N = 3 and Vq — 0. 
As was shown in Refs. P"7| , |T8fl , Eq. ( |1.1| ) describes, in fact, the transition of liquid helium-3 
from normal to superfluid Anderson-Morel phase; coupling constants go and Ao entering 
formulas of Refs. JTTyTEl ] are easily seen to be identical to uq and wq respectively. 



III. RG SERIES, RESUMMATION AND SYMMETRIES 



As was already mentioned, the static critical behavior of the model Eq. (|1 . 1| ) has been 
studied in three dimensions within one- and two-loop RG approximations P2"| . The taking 



into account of two-loop contributions to the /3-functions and critical exponents was found 
to change drastically the results of the lowest-order RG analysis. In particular, it alters 
the total number of fixed points and avoids the degeneracy of the 0(2A)-symmetric fixed 
point which is four-fold degenerate, for N = 2, within the parquette approximation. On the 
other hand, some of the numerical results obtained on the base of the resummed two-loop 
RG expansions do not obey some exact symmetry relations (see below). In such a situation 
three-loop calculations turn out to be very desirable. 

We calculate the /3-functions for the Hamiltonian Eq. (|1.1|) within a massive theory. The 
renormalized Green function Gn(p,m) and four-point vertex functions Un(pi,m,u,v,w), 
Vr(pj, to, u, v, w), Wr(pi, to, u, v, w) are normalized at zero momenta in a conventional way: 

G R \0,m) =m 2 , 
dG^{p,m) = 

Qpi p 2 =0 

Un(0,m,u,v,w) — mu , (3.1) 
Vr(0, to, u, v, w) = uiv , 
Wr(Q, to, u, v, w) = mw . 



One extra condition is imposed on the ip 2 insertion: 



T { n' 2 \p,q,m,u,v,w] 



p=q=0 



(3.2) 



The value of the one-loop vertex graph at zero external momenta including the factor (iV+4) 
is absorbed in u, v, w in order to make the coefficient for u 2 term in f3 u equal to unity. The 
/3-functions obtained are as follows: 



„ 2 4 , , , 2 , 2 

p u = u-u - 



N + 4 



(uv + uw + w 2 ) + 



27(N + A) 2 



(41 N + 95)u 3 + 200u 2 v + 200u 2 w 



+46uv 2 + {46 N + 216)uw 2 + 92uvw + lUvw 2 + {36 N + 72)u> J 

7TT- — ^ [(2.69789 N 2 + 54.94038 N + 99.82021)w 4 + (26.58751 N 

4(A^ + 4) 3 L V ' v 

+329.22770) [u\ + u 3 w) + (2.48756 N + 221.36225)(uV + 2u 2 vw) + (2.48756 N 2 
+ 155.55980 N + 470.42246)w 2 u> 2 + 50.50080(m; 3 + 3uv 2 w) + (34.28057 N 
+626.66599)mww 2 + (8.11011 N 2 + 125.31213 N + 311.16081)uiu 3 + 110.42034u 2 iu 2 
+ (1.95355 N + 216.93358)tw 3 + (-5.20190 N 2 - 0.62829 N + 95.22334)u? 4 



(3.3a) 



Pv = v 



NTl {3U + l V + 4w) + 27(N + 4) 2 



(23 N + 185)w 2 + 362m; + 52Auw 



' -2.50221 N 2 



+ 136v 2 + 380vw + (28 N + 180W 2 

v ' J 4(A^ + 4) 3 

+41.85390 N + 234.66699)w 3 + (-0.01437 N + 720.91540)m 2 v + (8.98498 JV 
+1015.38106) A + 579.33309m; 2 + 1575.28532mm; + (151.47423 N 
+780.92014)mw 2 + 157.45847*; 3 + 604.53412w 2 w + (6.49576 N + 753.08966)tw 2 
+ (-3.27046 N 2 + 13.63522 N + 284.67391)w 3 ^ 



(3.3b) 



A, 



w 



2 . N 2 

1 -—4 i3u + V + J w) + 27(N + 4y 



+92)uw + 28v 2 + 56vw + (36 - 8 N)w' 



(23 N + 185)w 2 + 200m; + (54 N 

' -2.50221 N' z + 41.85390 N 



4{N + 4) 3 



+234.66699)w 3 + (-9.01372 N + 426.44974)m 2 v + (2.99978 N 2 + 83.14193 N 
+230.13930)m 2 w + 162.71394m; 2 + (29.26715 N + 266.89358)uuiu + (5.75601 N 2 
+48.11146 N + 131.38337)mi> 2 + 25.29977v 3 + (1.15422 N + 73.59085> 2 u> 
+ (9.52258 N + 106.38551>«; 2 + (-1.31497 N 2 + 10.71074 N + 58.66955)u> 3 



(3.3c) 
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Such series are known to be divergent, at best asymptotic. They contain, however, 
rich and important physical information which may be extracted provided some procedure 
making them convergent is applied. The Borel transformation usually plays a role of this 
procedure. Here we are dealing with expansions of quantities depending on three variables 
u, v, and w. So, Borel transformation should be taken in the generalized form: 

oo 

f(u,v,w) = Y^c ijk u t v j w k = J e^F^tjVtjW^dt , (3.4) 

where the Borel transform expansion is as follows: 

^Hgf-^f^W (3.5) 

To calculate the integral entering Eq. (|3.4|) one should perform an analytical continuation 
of the Borel-transformed expansion. Although there are several different ways to do it, only 
two approaches proved to be efficient in the case of multi-variable RG series [f|||. The first 
one exploits the so called resolvent series [|47|: 



oo n n—l I m y n—l—m 

F(x, ^,A) = EA n EE Q ' m ' n "'- m 7 (3-6) 

n=0 1=0 m=0 U - 

which is actually a series in powers of A with coefficients A n being uniform polynomials of 
n-th order in u, v, and w. Pade approximants in A [L/M] is then used and the sum of the 
series is given by 

F(x,y,z) = [L/M]\ x=i (3.7) 

(see Appendix for detail). This approximation scheme possesses the remarkable property: 
for y = z = 0(orx = z = 0, OTX = y = 0) expression (|3.7f ) turns into conventional 
single-variable Pade approximants. Hence, all the results obtained for simpler, say, 0(n)- 
symmetric models hold good within this approach. 

Another way of the analytic continuation is realized through the construction of the 



Canterbury approximants invented by Chisholm ||48fl : 

K L M 

EE E A klm x k y l z m 

[K,L,M/R,P,Q} = k =° l =° p m = Q ° . (3.8) 

EEE B rpq x r yPzi 

r=0 p=0 q=0 
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It was found to be rather effective when applied to the impure Ising model [0,0, the cubic 
model f|||, and the mn-model 0. 

To determine which approximation scheme is the most adequate to our problem certain 
criterious should be formulated. We adopt the following ones: 

i) the resummation technique chosen should not lead to unphysical results, 

ii) new results should be consistent with the most accurate numerical estimates for 0(n)- 
symmetric and other simple models known up today, 

iii) new results should be self-consistent, i.e. numerical values of any critical exponent 
calculated by means of the resummation of different expansions, say, expansions for 7 and 
7 _1 , should be identical (as close as possible), 

iv) all (known) symmetries of the problem should be preserved by the approximation 
scheme employed. 

The last criterion is of prime importance in the case considered. The point is that the 
model Eq. ( |1.1| ) for N = 2 possesses specific symmetry properties. Indeed, if the field <p a 
undergoes the transformation 

fx -> <Pi , V?2 -> i<f2 , (3.9) 

quartic coupling constants are also transformed: 

u — > u , v—*v + 2w , w — > — w , (3.10) 

but the structure of the Hamiltonian itself remains unchanged ||2] . Just the same situation 
takes place in the case of another field transformation |n,l3"f: 



91 ~jr~ ' V2 ^^jt~ ' (3 - n) 

which does not affect the Hamiltonian structure resulting only in the following replacement 
of u, v, and w: 

v 

u — y u + v + 2w , v — > — 2w , w — > — — . (3.12) 

It is well known that RG functions of the problem are completely determined by the 
structure of the Hamiltonian: they do not depend on uq, Vq, and Wq which play a role of 
initial values of effective coupling constants when the RG flow of u, v, and w is searched. 
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Hence, RG equations should be invariant with respect to any transformation conserving the 
structure of the Hamiltonian ||49| ; Eqs. (|3.10| ) and ( 3.12j) , in particular, were shown to be 
such transformations. 

It means that for N = 2 /3 U , j3 v , and (3 W should obey some special symmetry relations 
which may be readily written down: 



(3 u (u, v, w) = [3 u (u,v + 2w, -w) , 

(3 v {u, v, w) + 2(3 w (u, v, w) = f3 v (u, v + 2w, —w) 

P w (u,v,w) — — P w (u,v + 2w, — w) . 

and 



(3.13) 



f3 u (u, v, w) + (3 v (u, v, w) + 2ft w (u, v, w) = 

V 

= p u (u + v + 2w, -2w, --) , 

V 

P v (u, v,w) = -2(3 w (u + v + 2w, -2w, --) , 

V 

2(3 w (u,v, w) = —(3 v (u + v + 2w, -2w, —-) . 



(3.14) 



One can see that expansions Eqs. (|3.3a| , |3.3b| , p.3c|) do really satisfy these relations. Moreover, 
due to this special symmetry, transformations Eqs. ( |3.10|) and ( |3.12|) can, at most, rearrange 
the fixed points of RG equations not affecting numerical values of their coordinates u c , v c , and 
w c themselves. It provides powerful tool for evaluation of the accuracy of the approximation 
scheme employed. 

To calculate critical exponents field-theoretical expansions for two of them are needed. 
We find 7 _1 and r] as power series in u, v, and w up to three-loop order. They are as follows: 
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iV + 4 



N + l 



-u + V + w 



+ 



(N + A)' 



N + l 



u 2 + 2{uv + uw + vw) + v 2 + Nw' 



0.2472701 



(N 2 + 5N + 4)w 3 + (6iV + 24)(u 2 v + u 2 w + vw 2 ) + 10(3uv 2 + 6uvw 

0.1925093 



(iV + 4) 3 

W + 3v 2 w) + (18N + 12)uw 2 + (2N 2 + 8)w 3 
18(JV+1) 



(N +2N + IK 



(N + 4) 3 

(u 2 v + u 2 w) + (2A^ + 10) {uv 2 + 2uvw) + (2N 2 + 2N + 8)uw 2 



+4{v 3 + 3v 2 w + Nw 3 ) + (4N + 8)vw' 



(3 
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27(A + 4) 2 
0.01234194 r 



(N + l)u 2 + 2(2uv + 2uw + v 2 + 2vw + Nw 2 ) 



+- r^- (N 2 + 5N + A)u 3 + (6N + 24)(u 2 v + u 2 w + vw 2 ) 

(iV + 4) 3 L v ' 

+10(3m; 2 + 6uvw + v 3 + 3v 2 w) + (18 N + 12)uw 2 + {2N 2 + 8)w 3 



(3. 



Since critical exponents are measurable (observable) quantities, the right-hand sides of 
Eqs. ( |3.15| ) and ( |3.16| ) should contain for N = 2 only those combinations of u, v, and 
w which are invariant under the transformations Eqs. ( |3.10|) and (|3.12 ). As may be seen, it 
is actually the case. 



IV. NUMERICAL RESULTS 



So, we perform the resummation of three-loop expansions Eqs. ( |3.3a| , |3.3b| , |3.3c| ) by 



means of generalized Pade-Borel technique with the approximant [3/1] being used for the 
analytic continuation of Borel transforms. Coordinates of the fixed points of RG equations 
thus obtained are found numerically for the most interesting cases N = 2 and N = 3. 
They are presented in Table | (N = 2) and Table [□] (N = 3) which contain also, for 
comparison, the fixed point coordinates obtained earlier from two-loop RG expansions 
resummed on the base of [2/1] Pade approximants. Three-loop contributions are seen to 
change appreciably locations of all non-trivial fixed points. 

Let us first discuss the numerical accuracy of the values found. The point 2 in Tables | 
and J| is actually 0(2A)-symmetric fixed point and its coordinates are to be compared with 
those obtained for 0(4)- and 0(6)-symmetric models with real fields (p a from resummed 
highest-order RG series available up today. Four-loop calculations in 3D have resulted in 
u c = 1.377 for n = 4 and u c = 1.338 for n = 6 ||50|| . These numbers differ from their 
three-loop counterparts presented in second columns of Tables | and [Tl] by no more than 
1%. 

Third columns of the Tables | and [TT] contain coordinates of the Bose {XY) fixed point. 
The most accurate estimate for v c obtained by the resummation of the six-loop 3D RG 
series is v c = 1.405 To compare this number with those presented in the Tables | 

and m however, we should make some rescaling of v c for N = 2 and N = 3. The point is 
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that the coefficient for v 2 in (3 V (Eq. ( |3.3b|) ) is equal to j^rg differing from unity for iV 7^ 1. 
Since the six-loop value of v c has been calculated in 0(2)-symmetric model, i.e. for N = 1, 
the numbers in third columns of Table | and |1| should be multiplyed, before comparison, 
by factors § and § respectively. It gives v c = 1.4032 (N = 2) and v c = 1.4033 (N = 3). 
Practical coincidence of these two values is very natural since they are actually coordinates 
of the same (Bose) fixed point while their closeness (~ 0.1%) to the six-loop value of v c 
provides an evidence of the high accuracy of the approximation scheme employed. Note, 
that the two-loop approximation leads to v c = 1.5583 which is more than 10% away from 
the "exact" value. 

Strong evidences of the high numerical accuracy of the approach elaborated may be ob- 
tained on the base of symmetry arguments. As was shown above, transformations Eqs. (|3.10|) 



and ( |3.12| ) can only rearrange the fixed points of RG equations ( |3.3a| , |3.3b| , ^.3c| ) for N = 2 



not affecting the values of u c , v c , and w c themselves. Indeed, this is precisely what occurs 
when one applies Eq. ( |3.10| ) to the content of Table [IJ the first four fixed points stay at their 



places while the points 5-8 undergo pair transpositions 5 — > 7, 7 — > 5, 6^8, 8^6. An- 
other transformation, Eq. ( |3.12| ), practically does not change the location of fixed points 1, 
2, 7, and 8 and causes pair transposition 3 — > 5, 5 ^ 3. The rest of fixed points, the 4-th and 
the 6-th ones, however, are converted one to another under Eq. ( |3.12| ) only within three-loop 
approximation. Corresponding two-loop results turn out to violate the symmetry relations 
induced by Eq. (|3.12|) . More precisely, the differences between the coordinates of the point 
4 and the transformed coordinates of the point 6 ("symmetry discrepancies") given by [2/1] 
Pade-Borel approximants are about 0.3, while within three-loop approximation they are of 
order of 0.01. 

So, the three-loop terms being taken into account enable one to obtain results which 
are much more accurate than those given by two-loop RG expansions. Moreover, it is 
seen that the field-theoretical RG approach in three dimensions combined with a general- 
ized Pade-Borel resummation technique does really provide a regular, rapidly converging 
approximation scheme powerful enough to treat complicated model with three quartic cou- 
pling constants. At the same time, the Chisholm-Borel resummation procedure is found to 
give poor results in this case. 

Let us discuss further the stability of the fixed points and the structure of the RG 
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flow diagrams. All fixed points of the RG equations are unstable in the three-dimensional 
parameter space (it, v, w). The 4-th and the 6-th ones, however, are stable within the planes 
(u, v) and (it, w) respectively. The existence of such points is important since it implies the 
possibility of continuous phase transitions in numerous physical systems described by the 
model Eq. (|1 . 1| ) with N = 2 and Vq — or iy = 0. RG flows for N = 2 within the planes 
(u, v) and (u, w) and for N = 3 within the plane (it, v) are shown in Fig. [TJ(b) and Fig. |2](b, 
c). One can see from these figures that there is not a fixed point stable within the plane 
(u, w) for N = 3 while for N = 2 such a point exists. Hence, the topology of the flow 
diagram should change when N varies. It is interesting, therefore, to study the structure of 
our RG flows for arbitrary N. 

The detailed numerical analysis of three-loop RG equations obtained shows that only two 
diverse u — v flow pictures occur for 1 < iV < oo while the RG flow within the plane (u, w) 
may proceed in four different ways. All possible scenarios are depicted in Fig. [I] and Fig. |2| 
(see hardcopy: Phys. Rev. B 49, 15901 (1994)). The critical (marginal) dimensionality of 
the order parameter N c which separates from each other two regimes of RG flows for w — 
is found to be: 

N c = 1.47 ±0.01 . (4.1) 

Since this number is less than two, in all physically interesting cases, i.e. for N > 2, the 
0(2A r )-symmetric fixed point turns out to be unstable. So, the system should undergo 
either a continuous phase transition demonstrating an anisotropic (v c ^ 0) critical behavior 
or a fluctuation induced first-order phase transition. When iV — > oo the anisotropic stable 
fixed point in the plane (it, v) is going to the 0(2)-symmetric one which becomes degenerate 
in this spherical-model limit. 

The behavior of our model in the plane (u, w) is more rich. It is characterized by three 
marginal values of the order parameter dimensionality: iV cl , and N c ^. Calculated on 
the base of resummed three-loop RG series Eqs. ( |3.3a| , |3.3b| , |3.3c| ) they are as follows: 



N cl = 1.45 ±0.01 , 

N c2 = 2.03 ± 0.01 , (4.2) 
N c3 = 3.91 ± 0.01 . 
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For iV < N c i (Fig. 0(a)) three non-trivial fixed points exist in the plane (u,w) with the 
0(2iV)-symmetric point being stable. When N exceeds iV cl this "Heisenberg" fixed point 
loses its stability but the other, anisotropic fixed point with w c < acquires it (Fig. @(b)). 
In this domain which includes the important case N = 2, our system demonstrates an 
anisotropic scaling behavior or discontinuous phase transitions. With increasing N the stable 
fixed point in Fig. 0(b) is moving downward and "annihilates" with the anisotropic saddle 
fixed point when N approaches N C 2. There is only one non-trivial fixed point in the domain 
N C 2 < N < N c3 including N = 3 (Fig. 0(c)); it is 0(2iV)-symmetric and unstable. So, only 
first-order phase transitions are possible, in principle, in this case provided wq ^ 0. At last, 
when N increases further and crosses over the value N C 3 the creation of two new anisotropic 
fixed points in the u — w flow diagram takes place (Fig. @(d)). One of them is stable and 
describes some anisotropic critical behavior with w c > 0. This fixed point is known as a 



"chiral" point |L3[ and corresponding "chiral" phase transition has been extensively studied 
during the last years. As follows from our estimates (Eq. ( (4.2[ )), this point does really exist 
and governs the scaling behavior of physical systems with N > 4. For N = 2 and N = 3 
the chiral critical behavior does not actually realized. 

Let's discuss numerical estimates Eqs. (4J.) and (|4.2j ) in more detail. The value of N c2 



turns out to be very close to iV = 2 which is of prime physical importance. Can higher-order 
contributions to the /3-functions being taken into account change iV c2 , invert the inequality 
N C 2 > 2, and, hence, alter the structure of the u — w flow diagram for N = 21 No, they 
can not. The point is that the structures of the RG flows in the planes (u, v) and (u, w) are 
related to each other for N = 2 by the symmetry relations discussed earlier. In particular, as 
may be seen from Eq. ( |3.12|) the total number of fixed points in each of these flow diagrams 



should be just the same. Since the plane (u, v) definitely contains four fixed points (iV c lies 
far below the value of interest N = 2) the plane (u, w) for N = 2 should possess four fixed 
points too. Moreover, since, for N = 2, the stable fixed point has v c > its counterpart in 
the plane (u, w) should possess w c < (see Eq. fl3.12|) ). It means that inevitably N c2 > 2 in 



the exact theory. The estimate Eq. ( |4.2|) is in accord with this inequality. 

Another point to be discussed is the near coincidence of the calculated values of N c 
and N c i. It is not occasional. Indeed, N c and iV cl are both the values of iV for which 
the 0(2iV)-symmetric fixed point becomes degenerate and critical exponents describing its 
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stability change a sign. These exponents are completely determined by derivatives ^ and 
taken at the "Heisenberg" fixed point since ^ and ^r- at this point vanish. One can 

aw ° r aw av r 

see from Eqs. ( |3.3b| ) and ( |3.3cj ) , however, that ^ = 4^ along the whole line v — w — 
up to the highest calculated order. So, when iV varies the "Heisenberg" fixed point should 
lose its stability in planes (u, v) and (u, w) simultaneously, i.e. N c and iV cl should be equal 
to each other. The small difference between calculated values of N c and iV cl reflects a finite 
accuracy of our approximation scheme which is seen to be of order of 0.01. 

Having calculated fixed point coordinates we can find critical exponents for our model. 
To obtain accurate estimates for 7 the expansion Eq. ( |3.15| ) is resummed by means of the 
generalized Pade-Borel procedure described above while the values of 77 are found by direct 
substitution of fixed point coordinates into Eq. ( |3.16| ) since this very short series with very 
small and positive three-loop term doesn't need in resummation. The results obtained for 
N = 2 and N = 3 are presented in Table [Hj and Table [iy| respectively which contain also 
the values of 7 and rj calculated earlier ]42[ within two-loop approximation. 

Three-loop contributions are seen to change critical exponents values only slightly. For 
N = 2 critical exponents calculated in the fixed points 3, 5, and 7 turn out to be almost 
identical, and so is true for the fixed points 4, 6, and 8. In the exact theory each of these 
two sets of fixed points indeed should possess identical critical exponents since fixed points 
belonging to the same set are related to each other by symmetry relations Eqs. ( |3.10| ) and 
( ^.12| ), i.e. they are actually the same fixed point. So, differences between the values of 7 and 
r) calculated in such fixed points may be considered as a measure of numerical accuracy of our 
approximation. It is seen to be better than 0.001. On the other hand, the difference between 
the values of 7 calculated in the fixed points 4 and 6 (or 8) within two-loop approximation 
exceeds 0.05. Hence, the taking into account of three-loop terms improves the situation 
essentially. 

It is worthy also to compare critical exponents found in the Bose and "Heisenberg" fixed 
points with their counterparts determined from six-loop |l|J^] and four-loop RG expan- 
sions for a simple 0(n)-symmetric model. The most accurate estimate for the susceptibility 
exponent of the 3D XY model is 7 = 1.315 [0]||. Corresponding values in Tables |T| and 
P"V| (third columns) differ from it by 0.005. Four-loop RG calculations for n = 2N = 4 and 
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n = 2N = 6 give 7 = 1.441 and 7 = 1.541 respectively fl50| . Differences between these num- 
bers and their three-loop twins presented in Tables |T| and fV] (second columns) are about 
0.02. So, we arrive to the conclusion: Pade-Borel resummed 3D three-loop RG expansions 
provide an accuracy of order of 0.01 for all calculated quantities. This accuracy is sufficient 
for making definite and reliable theoretical prediction for physical systems described by the 
model Eq. ( |Llp . It will be done in the following Section. 

Now let's return back to the calculation of critical exponents. The rest of them may be 
found by making use of well-known scaling relations. We present here numerical values of 
exponents u, a, and (3 for the fixed points which are stable within corresponding parameter 
subspaces since only these numbers may be related to experiments. So, for equivalent fixed 



points 4, 6, and 8 from Table LI I 



v = 0.677 , a = -0.030 , (3 = 0.347 (4.3) 



while the point 4 in Table [TV] is characterized by v = 0.673, a = —0.020, and (3 = 0.345. 



V. APPLICATION TO PHYSICAL SYSTEMS AND DISCUSSION 

All fixed points of our RG equations were found to be unstable within the three- 
dimensional parameter space (u, v, w) for N = 2 and N = 3. It means that only dis- 
continuous, first-order phase transitions should occur in physical systems with non-zero 
initial values of v and w. Such systems are represented by cubic and tetragonal unconven- 
tional superconductors and superconductors with composite s — d order parameters. On 
the other hand, fluctuation-driven first-order phase transitions are known to be extremely 
weak. So, the absence, within the experimental accuracy, of discontinuous superconducting 
transitions in relevant heavy-fermion and high-T c compounds does not actually contradict 
to the above conclusion. 

In hexagonal ci-wave superconductors described by the model Eq. ( |1 . 1| ) with N = 2 
and v = second-order phase transitions remain possible under strong superconducting 
fluctuations since there is a stable fixed point within the plane (u, w) which possesses a 
sizable domain of attraction. Corresponding values of critical exponents (column 4 (6, 8) 



in Table III and Eq. (13)) turn out to be close enough to those of the 3D XY model. 
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So, it is actually impossible to distinguish between the BCS s-wave pairing and the non- 
trivial one studying experimentally the scaling behavior of superconductors. On the other 
hand, anisotropic gradient terms omitted in the Hamiltonian (|1 . 1| ) can themselves change, 
in course of fluctuation renormalization, the order of phase transition and the structure 
of phase diagram, as they do in crystals undergoing structural (ferroelectric) and magnetic 
phase transitions ]5T| , |52| . This will obviously result in a non-universal behavior of hexagonal 
<i-wave superconductors in the critical region. 

In liquid helium-3, where N = 3 and v o = and, therefore, RG equations have no stable 
fixed points, fluctuations should always force the superfluid phase transition to be first- 
order. Corresponding discontinuities of thermodynamic quantities at the transition point, 
however, would hardly be observed experimentally because of the narrowness of the critical 
region in this Fermi-liquid (see, e.g. Refs. f° r numerical estimates). 

Only first-order phase transitions should emerge also in FFJJA at T = 0, in spite of the 
existence of stable fixed points in planes (u, v) and (u, w) for N = 2. Indeed, RG trajectories 
starting from physical initial points, i.e. from those having v o < and Wq > for triangular 
and square FFJJA respectively |TT| , |4"5|| , can not achieve the stable fixed points as is cleary 
seen from Fig. |I](b) and Fig. ||(b). So, these systems will demonstrate non-universal critical 
behavior. 

A mode of the critical behavior of STA and HM described by the Hamiltonian ( |1.1|) with 



v o = depends on the dimensionality of the order parameter. In materials with Heisenberg 
spins, i.e. for N = 3, only (weak) first-order phase transitions should occur. In easy-plane 
crystals with Xy-like spins continuous transitions are also possible with critical exponents 
presented in column 4 (6, 8) of Table [H| and Eq. (|4.3| ) which are practically undistinquishable 
from those of 3D XY model. These exponents, however, govern transitions into somewhat 
trivial phases: simple ferromagnetic or antiferromagnetic in STA and a sinusoidal (linearly- 
polarized) in HM, since relevant stable fixed point possesses w c < 0. Much more interesting 
ordering with frustration in STA and helical one in HM are described by Eq. ( |1 . 1|) with 
wq > 0. They may be realized only via first-order phase transitions, as is clearly seen from 
Fig. |(b). 

We did not find any traces of chiral second-order phase transitions and corresponding 
new classes of universality for N = 2 and N = 3, i.e. for STA and HM with Heisenberg or 
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XY-\ike spins. This result is in contradiction with conjectures and conclusions made on the 
base of the lower-order e-expansion analysis |l|f5|. Such conclusions, however, can not be 
referred to as reliable since the method mentioned provides rather low numerical accuracy in 
three dimensions. To illustrate this point and to clear up the situation let's discuss two-loop 
e-expansions (highest-order now available) for marginal order-parameter dimensionalities 
N&, N c2 , and N c3 . They are as follows 



iV cl = 2-e , (5.1a) 
N c2 = 2.20 - 0.57e , (5.1b) 

7V c3 = 21.8 - 23.4e . (5.1c) 

When e — > 1 N C 3 becomes less than 2 and chiral fixed point seems to exist for N = 2 and 
N = 3. In this limit, however, N c2 becomes less than 2 too what is in obvious contradiction 
with an inequality N c2 > 2 proven above. Moreover, another inequality N c2 < N c3 valid for 
e < 1 turns out to be broken at e = 1 as well. 

Is it possible to make e-expansion predictions more accurate for e = 1? Yes, of course. 
Higher-order (four- and five-loop) e-expansions are known to give rather good numerical 
results for 3D 0(n) -symmetric model at e = 1 provided some Borel-like resummation 
procedure is applied [[51H5"6|1- Unfortunately, we have no long enough e-expansions for our 
model up today |57fl . So, we try to "sum up" expansions ( |5.1a| , |5.1b| , |5.1cj ) constructing 



simple Pade approximants: 



iV cl « , (5.2a) 

l + 0.5e ' V ' 



2 2D 

N c2 « , (5.2b) 

l + 0.26e v ; 



"c3 



21.8 

(5.2c) 



1 + 1.07e 



For e = 1 these formulas give iV cl = 1.33, N c2 = 1.75, and N c3 = 10.5. The first number 
is much closer to our estimate N c i = 1.45 (Eq. Q4.2D ) than the value N c \ = 1 given by 
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Eq. ( |5.1a| ). The second one is also closer to the three-loop 3D estimate N c2 = 2.03 than the 
naive value N c2 = 1.63, but both violate the inequality N c2 > 2. The third number exceeds 
enormously the estimate N c3 = 3.91 which turns out to lie between this number and the 
naive estimate N C 3 = —1.6. So, we see that primitive resummation of very short expansions 
( f>. la , [5.1b , |5.1c|) results in somewhat improved numerical estimates for iV cl and N c2 while 
being used for evaluation of N c3 it demonstrates that lower-order e-expansions are useless 
in this case. Hence, lower-order calculations in 4 — e dimensions can not be considered as 
a tool for answering the question whether 3D physical systems with N = 2 and N = 3 
undergo chiral phase transitions or not. 

Monte-Carlo simulations UT2"| , |2"2"|| would also hardly be referred to as evidence of chiral 
critical behavior of STA and HM with Heisenberg or Xy-like spins. The point is that 
unusual values of critical exponents given by such calculations turn out to be close to tricrit- 



ical ones. That's why it was suggested |23| that tricritical behavior or tricritical-to-critical 
crossover are really seen in these computer experiments as well as in most of physical ex- 
periments performed on several helimagnets (Tb, Dy, Ho) and STA (CsMnBr 3 , CSVCI3 and 
others). We completely agree with what is argued on this topic in Ref. where the reader 
can find also an overview and analysis of relevant experimental data. 



VI. CONCLUSIONS 



The critical behavior of the model describing phase transitions in superconducting and 
magnetic systems with complex iV-vector order parameter as well as in superfluid helium- 
3 has been studied within the RG approach in three dimensions. RG /3-functions and 
critical exponents have been calculated as series in powers of renormalized quartic coupling 
constants u, v, and w up to three-loop order. The series obtained have been resummed 
by means of the generalized Pade-Borel technique and fixed points coordinates, critical 
exponents values and a structure of RG flows have been determined for N = 2 and N = 3. 
Marginal values of the order parameter dimensionality at which the topology of RG flows 
in planes (u, v) and (u, w) changes have been also found. Several criterions have been used 
to estimate the accuracy of numerical results obtained which had turned out to be about 
0.01, an order of magnitude better than that given by resummed two-loop RG expansions. 
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So, field-theoretical RG approach in three dimensions combined with a proper resummation 
technique provides a regular, rapidly converging approximation scheme powerful enough to 
treat complicated model with three quartic coupling constants. 

Relevant physical systems have been shown to undergo, in most cases, fluctuation induced 
first-order phase transitions. Second-order transitions have been found to occur only in 
hexagonal c?-wave superconductors and in planar magnets (into somewhat trivial phases: 
linearly-polarized or unfrustrated) . Corresponding critical exponents have turned out to 
differ from those of the 3D XY model by no more than 0.02 - 0.03, i.e. the underlying 
critical behavior would hardly be thought as experimentally distinguishable from the Bose 
one. RG equations obtained have been shown to possess the chiral fixed point but only for 
N > A. It means that STA and HM with Heisenberg and AF-like spins would not really 
demonstrate the chiral critical behavior with unusual critical exponents approaching helical 
or frustrated antiferromagnetic states via first-order phase transitions. 
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In this Appendix, some details of the resummation procedure employed are described. 
As was shown in Sec. pi] , resolvent series 

oo 
n=0 

(1) 

n n—l 

A n = ^2^2 ' m ' n ~ l ~ m x l y m z n ~ l ~ m 

1=0 m=0 n ' 
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for Borel transforms of the original multi-variable RG expansions may be constructed to 
generate Pade approximants [L/M] in parameter A. These approximants are defined in a 
conventional way: 

{L/M] = £i ■ (2) 

where -Pl(A) and Qm(A) are polynomials of degrees L and M respectively with coefficients 
depending on x, y, and z, which may be determined from the following relations: 

Q M (X)F(x,y,z; A) - P L (X) = 0(A i+M+1 ) , 

(3) 

Qm(0) = 1 . 

Approximate expressions for /3-functions and critical exponents are then obtained by the 
replacement of variables x = ut, y = vt, and z = wt in Pade approximants and by evaluation 
of the integral 

oo 

Je-'iL/M^^dt (4) 
o 

(Borel transformation). 

With three-loop expansions in hand, we can use two different approximants [3/1] and 
[2/2] obeying the condition L > M. The former was shown (Sec. |TV| ) to provide rather good 
numerical results for all cases considered. Moreover, an employment of this approximant 
kept us away from the well-known problem of poles which often arises when approximants 
with higher-order denominators are used. That's why we have chosen Pade approximant 
[3/1] for our analysis. When expressed in terms of renormalized coupling constants u, v, 
and w and the variable t it is as follows: 

rq/ii _ Q o + ait + a,2t 2 + a^t 3 . . 

[ 1 1 ~ l + M ' 

where a , . . . ,a 3 and b% are known functions of u, v, and w. If series to be resummed are 
those for /3-functions the coefficient a in Eq. (]5|) turns out to vanish and the integral (JD 
reads 



/ 



te _ t a 1 + a 2 t + a 3 t* dt 



l + M 
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Evaluating this integral we get the final expression (the "sum" of the series) for the function 
of interest: 



f(u, v, w) = (ai + a 2 + 2a 3 )6 - (a 2 + a 3 — a 3 b)b 2 
+ (ai - a 2 6 + a 3 b 2 )b 2 e b Ei(-b) , 



(7) 



where Ei(x) is the exponential integral [5B| and b = b x l . This is precisely the formula 
which was used for resummation of the three-loop RG expansions Eqs. ( p. 3a , p. 3b , |3.3c| ) 



and for determination of the fixed points. The approximate expression for 7 
quite similar and not presented here. 
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FIGURES 

FIG. 1. RG flows in the plane (u, v) for N < N c and N > N c where N c = 1.47 ± 0.01. Shaded 
areas represent the regions of instability of the Hamiltonian ( |1.1|) (see hardcopy: Phys. Rev. B 
49, 15901 (1994)). 

FIG. 2. Four possible scenarios of RG flow in the plane (u,w). Marginal values of the order 
parameter dimensionality N c i, N C 2, and N c % are given by Eq. ( |4.2| ). Shaded areas are the regions 
of instability of the Hamiltonian flLl| ) (see hardcopy: Phys. Rev. B 49, 15901 (1994)). 
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TABLES 



TABLE I. 


Coordinates of the fixed points of 


RG equations 


for N = 


2 obtained 


within 


three-loop (approximant [3/1]) and two-loop (approximant [2/1]) approximat 


ons. 






i 
i 


O 
Z 


o 


1 


c 



D 


7 


8 


Ur f3/ll 


0.0 


1.3671 


0.0 


0.1872 


1.6833 


1.6787 


1.6832 


1.6789 


[2/1]* 


0.0 


1.4863 


0.0 


0.0340 


1.8699 


1.8334 


1.8699 


1.8334 


v c [3/1] 


0.0 


0.0 


1.6838 


1.4914 


0.0 


0.0 


-1.6800 


-1.4950 


[2/1]* 


0.0 


0.0 


1.8699 


1.8334 


0.0 


0.0 


-1.8699 


-1.3591 


w c [3/1] 


0.0 


0.0 


0.0 


0.0 


-0.8416 


-0.7477 


0.8400 


0.7480 


[2/1]* 


0.0 


0.0 


0.0 


0.0 


-0.9350 


-0.6796 




u.u i yo 


* Quoted from Ref. [42 
















TABLE II. 


Coordinates of the fixed points of RG equations for N = 


3 obtained 


within 


three-loop (approximant [3/1]) and two-loop (approximant [2/1]) approximat 


ons. 










1 




2 


3 




4 


u c 


[3/1] 




0.0 


1.3310 


0.0 




0.0780 




[2/1]* 




0.0 


1.4262 


0.0 




0.0097 


v c 


[3/1] 




0.0 




0.0 


1.9646 




1.8845 




[2/1]* 




0.0 




0.0 


2.1816 




2.1713 


w c 


[3/1] 




0.0 




0.0 


0.0 




0.0 




[2/1]* 




0.0 




0.0 


0.0 




0.0 



* Quoted from Ref. §2| 
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TABLE III. Critical exponents 7 and rj for N = 2 calculated within three-loop (approximant 
[3/1]) and two-loop (approximant [2/1]) approximations. 





1 


2 


3 


4 


5 


6 


7 


8 


7 [3/1] 


1 


1.4260 


1.3099 


1.3360 


1.3098 


1.3355 


1.3102 


1.3357 


[2/1]* 


1 


1.4347 


1.3218 


1.3259 


1.3218 


1.3799 


1.3218 


1.3799 


V [3/1] 





0.0257 


0.0261 


0.0261 


0.0260 


0.0261 


0.0260 


0.0261 


[2/1]* 





0.0273 


0.0288 


0.0287 


0.0288 


0.0286 


0.0288 


0.0286 


* Quoted from Ref. [42 


1 














TABLE IV. Critical exponents 7 and r] for A = 3 calculated within three-loop (approximant 


[3/1]) and two-loop (approximant [2/1]) approximations. 
















1 


2 




3 




4 


7 


[3/1] 




1 


1.5164 




1.3099 




1.3291 




[2/1]* 




1 


1.5217 




1.3218 




1.3220 


V 


[3/1] 







0.0238 




0.0261 




0.0261 




[2/1]* 







0.0246 




0.0288 




0.0286 



* Quoted from Ref. || 
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